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Abstract. A local mesh refinement procedure for the CE/SE method which does 
not use an iterative procedure in the treatments of grid-to-grid communications is de- 
scribed. It is shown that a refinement ratio in the order of 10 can be applied successfully 
across a single coarse grid/fine grid interface. 


1 Introduction 

To accurately resolve flow properties in high- gradient regions without applying 
a very fine and costly uniform grid over the entire computational domain, simu- 
lation of practical unsteady CFD problems may require implementation of local 
mesh refinement (LMR) in both spatial and temporal directions. In this paper, 
we report a new LMR procedure for the CE/SE method [1-3]. 

It is a must for a LMR procedure to enforce space-time flux conservation 
across a coarse grid/fine grid interface so that numerical solutions wdll not be 
contaminated by spurious reflections originating there [4]. Fortunately, this re- 
quirement can be met easily in the CE/SE setting. In particular, an iterative 
procedure is not needed in the treatments of grid-to-grid communications. Fur- 
thermore, it will be shown that a refinement ratio in the order of 10 can be 
applied successfully across a single interface. 


2 Preliminaries 


To simplify the description of the LMR procedure, consider the PDE 

du/dt + adujdx = 0 (1) 


and its conservation form 


/ 

JS(V) 


h - ds = 0 


( 2 ) 


where (i) a > 0 is a constant; (ii) S(V) is the boundary of any region V in 
space-time; (iii) h = (au, u); and (iv) ds = dcrti with dcr and n, respectively, 
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being the area and the outward unit normal of a surface element on ^(V). Note 
that h * ds is the flux of h leaving V through the surface element ds. 

Let J?i denote the set of all mesh points (j, n) (dots in Fig. 1) with j + n 
being odd integers. For each (j,n) € ft i, let the solution element SE(j,n) be the 
interior of the space-time region bounded by a dashed curve depicted in Fig. 2). 
For any (x, t) € SE(j, n), let u(x,t) and h(x, t), respectively, be approximated 
by u*(x,t]j,n) and h where [1] 

u*(x,t;j,n) = uj + (u x )j(x-Xj) + - t n ) (3) 

h*(x,t;j,n) = (au*(x,t;j,n), u*(x,t;j,n)) (4) 

It is required that u = u*(x,t\ j,n) satisfy Eq. (1). Thus 

(u t )] = -a(u x )” (5) 

Hereafter, qj 1 denotes the column matrix formed by u 7 - and (u x )", the indepen- 
dent marching variables at (j, n) € f2 \ . 

Let space- time be divided into conservation elements (CEs). Two CEs, i.e., 
CE_(j, n) and CE +(j,n) (Figs. 3-4) are assigned to each (j,n) € ft\. Because 
the bomidary of CE_ (j,n) (CE + (j,n)) lies within the union of SE(j,n) and 
SE(j - l,n — 1) (SE(j + l,n - 1)), q* 1 can be solved for in terms of q”±i by 
substituting V = CE±(j, n) and h = h* in Eq. (2). The resulting scheme is the 
a scheme which has the property that the flux of h* is conserved over the union 
of any combination of CEs [1], 

For the a-e scheme [1], which reduces to the a scheme if e — 0, again q 72 is 
a function of However, the flux of h* generally is not conserved over the 

CEs referred to above. Rather, for each (j, n) € J?i, it is conserved over CE(j, n) 
(the union of CE±(jf, n) — see Fig. 5). To proceed further, in the following, the 
concept of modified flux will be introduced. 

B y de finiti on, t he modified flux at each of the line segments AB, BC \ CD 

DE , EF , and FA, which form the boundary of CE(j,n), is the corresponding 
flux of h *. On the other hand, the modified flux at XD, i.e., the interface that 
divides CE ±(j, n), is defined such that the total modified flux leaving the bound- 
ary of CE_ (j, n) vanishes. Note that the modified flux at AD, except in the case 
of the a scheme, is not equal to the corresponding flux of h*. Because (i) the 
total modified flux leaving the boundary of CE( j,n) vanishes; and (ii) the modi- 
fied flux leaving CE_(j, n) through AD is the negative of that leaving CE+(j, n) 
through AD, one concludes that the total modified flux leaving the boundary 
of CE +(j,n) also vanishes. In other words, the modified flux is conserved over 
CE ±(j,n) at each (j, n) € ft\. 

The a-e scheme obviously can also be defined over J? 2 , the set of (j, n) (not 
marked in Fig. 1) with j + n being even integers. Moreover, two schemes defined 

over J?i and i ?2 can be fused into a “dual 1 ’ scheme which is defined over ft = f 
ft\ U i? 2 . Note that, in practical applications, a solution of a dual scheme may 
not be decoupled into two independent solutions. 
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Let (j, n) € l?i, then (j + l,n) € ft 2- By definition, CE+(j,?i) coincides with 
CE_(j -f l,n). However, tw r o different flux conservation conditions are defined 
over their co-boundary. In fact, two different modified fluxes tied to, respectively, 
these two coincident CEs, are defined at any boundary line segment. Obviously, 
the same complication arises for any CE. As a result, in case that the solution 
of a dual scheme becomes coupled, the concept of global flux conservation is 
not well defined in terms of these individual modified fluxes. However, if the 
generalized flux at any line segment that forms the boundary of a CE is defined 
as the ax r erage of the two local modified fluxes, then the generalized flux is 
conserved over the union of any combination of CEs. The proof follows from (i) 
the generalized flux is conserved over each CE^a fact that can be established 
by summing over the two local modified flux conservation conditions tied to this 
CE; and (ii) the generalized flux leaving one CE through an interface is the 
negative of that leaving the neighboring CE through the same interface. 


3 Time Marching 


As an example, a step-by-step description of the marching from the Oth to 1st 
time level (Fig. 6) will be described. Note that (i) the simple grid structure used 
here does not limit the generality of the following discussions; and (ii) two q s 
and thus two symbols (e.g., A\ and B{) are assigned to any common mesh point 
of a coarse grid and a fine grid. 

I. Because a > 0, q(A\) is specified as a part of the left boundary condition 
while q{Cz) is extrapolated from q(C§ ), be., q(C^) — q(C £)■ 

II. q{A\) is determined in terms of q( A ® ) and <7(A°) using the dual a-e scheme 
with e — 0.5 [1], Similarly, qlB^ 2 ) (^(Cj)) is determined in terms of q{B ® ) and 

3) (q( c \ ) and <?(C?)). 

III. q(B y 2 ) is determined in terms of <7(^2), <7(^3) and qiB®) assuming 


C (B, 1/2 ) = [u{B\ /2 ) - (u(A° 2 ) + {M/2)u t {A° 2 ))\ I AX 


F(A%B{ r *) + F(B[ /z B 2 /j ) + F{B%B i y) + F(B%B?) = 0 


73 


>1/2 r.1/3 


0 D 1 


% 


?0 \ 


( 6 ) 

(7) 


Here (i) u^A^) = (see Eq. (5)); and (ii) F(A^Bl^), as an example, 

denotes the flux of h* leaving the rectangle A^B^ 2 B^ 2 B$ through A^B^ 2 as 
evaluated using <7(^3). Similarly, q(B^ 2 ) is determined in terms of q(B §), qr(Cf) 
and qiC®)- <7(^2) then determined in terms of q(B 1 ^ 2 ) and q(B^ 2 ) using the 
special dual a-e scheme referred to above. 

IV. <7(i?J) is determined in terms of q{B^ 2 ), <7(^3) and <7(^2) assuming 


u x (Bl) = [u{B{) - (u(A%) + Atut(A°))] /ax 
F{A % , B[ /2 Bl) + F(B\B\) + F(B 2 /2 bI) + F(B^ / 3 B l y) = 0 


( 8 ) 

(9) 
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Here F(A®, B l ^ 2 B\) denotes the flux of h* at B\^ 2 B\ evaluated using q(A%). 
Similarly, q(Bl) is determined in terms of q(B ^ 2 ), q{Ci ) and q(C 2 ). 

V. q{A\) is determined in terms of q(Bf), q{B\ ^ 2 ), q(B \j 2 ) and q(A l§) as- 
suming 

= [(u(£ 2 1/2 ) + {At/2)u t {Bl /2 j) - u(4)] /(ax/ 2) (10) 

F(B°B^) + FiBY^Bl) + F(A lA*) + F(Ap|) + F(AA) = 0 ( 11 ) 

Similarly, q(Cl) is determined in terms of ^(^3), q(B^ 2 ), q(B^ 2 ) and q(C 2 ). 

Note that, at each line segment on a coarse grid/fine grid interface , two 
fluxes of h* are evaluated using a c.oarse-grid q and a fine-grid q, respectively . 
Moreover , each of these fluxes is involved in a coarse-grid conservation condition 
and a fine-grid conservation condition . As a result , for the current marching 
procedure , it can be shown that the generalized flux is conserved over the union 
of any combination of interior CEs. 

4 Mesh Reconstruction 

For time-dependent problems, spatial mesh reconstruction may be implemented 
at any designated time level such that high- gradient regions are always covered 
by fine grids. In the following, a step-by-step description of a reconstruction 
procedure that preserves generalized-flux conservation will be explained using 
Fig. 7 (Note: C3/2 is the midpoint of C1C2). 

I. No reconstruction is required for the mesh points to the left of A 2 or to 
the right of C2 . Thus q{A[) = q{A\) and q(C ' 2 ) = q(C$). 

II. q{A 2 ) is determined in term of q(A2), q(B\) and q(B 2) assuming 

F(jfi%) = F(ZZ 4 ?), F(A^ 3 ) = F(B^ 2 ) + F(B^) (12) 

III. q(B 2 ) is determined in terms of q(C\) and q(C%) assuming 

F(B&) = F(C 2 ,V^K), F(B$ 3 ) = F(C u C^&) (13) 

IV. <7(43) is determined in terms of <7 (#2)1 <7 (£3) and q(B 2 ) assuming 

u x (A' 3 ) = (u(B' 2 )-u(A' 3 ))/(ax/ 2), F{A£%) = F(B&) + F(Brf[) (14) 

V. q(B[) is determined in terms of q(C\) and q{A' 2 ) assuming 

u x (B' 1 ) = (u(B[)-u(A' 2 ))/ax, F(B[P 2 ) = F(U/C^ 2 ) (15) 

VI. <7(^3) is determined in terms of q(C2) and q{C 2 ) assuming 

u x (B' z ) = (u(C' 2 )-u(B' 3 ))/ax, F(B$ 2 ) = F(U&£) (16) 

VII. q(C[) is determined in terms of q{C2) and q{B 2 ) assuming 

u x (C[) = (u(C' 1 )-u(B' 2 ))/(ax/ 2), F(C[C?) = F(C^) (17) 
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5 Numerical Results 


Assuming a = 1, ax = 0.01 and At = 0.009, an exact solution to Eq. (1), i.e., 

u = (10/y/n) exp [— 100(x — t — 0.3) 2 ] (18) 

is solved numerically using the current LMR procedure with a grid refinement 
ratio R — 10 in both space and time. The initial condition (£ = 0) and mesh 
distribution are shown in Fig. 8 while the exact solution (a solid line) and nu- 
merical solution (dots) along with the mesh distribution at t = 0.9 are depicted 
in Fig. 9 (Note: due to a consideration related to generalized flux, here each dot 
represents a midpoint solution, i.e., the average mesh value of two neighboring 
mesh points). It is seen that the numerical solution match very smoothly across 
the two coarse grid/fine grid interfaces. 

Assuming ax = 4 x 10 -3 , At = 1.5 x 10 -3 and R = 4, the Sod’s shock tube 
problem [1] is also solved using a similar Euler LMR procedure. The results 
at t = 0.195 are shown in Fig. 10(a) while a close-up view covering only the 
neighborhoods of grid interfaces is depicted in Fig. 10(b). Note that a simulation 
using R = 64 also has been carried out successfully. 

6 Conclusions 

A simple CE/SE LMR procedure has been described. Its robustness is demon- 
strated by numerical examples with grid refinement in the order of 10. 
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Fig. 8. Initial condition (Eq. (18)). 
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